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Abstract 

The divide and color model on a graph G arises by first deleting each edge of G with 
probability 1 — p independently of each other, then coloring the resulting connected compo- 
£Nl ' nents (i.e., every vertex in the component) black or white with respective probabilities r and 

1 — r, independently for different components. Viewing it as a (dependent) site percolation 
D ' model, one can define the critical point r™ (p). 

In this paper, we first give upper and lower bounds for r^(p) for general G via a stochastic 
' comparison with Bernoulli percolation, and discuss (non-)monotonicity and (non-) continuity 

properties of r™ (p) in p. Then we focus on the case G = Z 2 and prove continuity of r J (p) as 
a function of p in the interval [0, 1/2), and examine the asymptotic behavior of the critical 
value function as p tends to its critical value. 

Oh 
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Introduction and main results 

> 

The divide and color (DaC) model is a natural dependent site percolation model introduced by 
Haggstrom in [5]. It has been studied directly in [5JEJ HIE]) and as a member of a more general 
family of models in HQ1 HI El E| • 

This model is defined on a multigraph G = (V, £), where £ is a multiset (i.e., it may 
contain an element more than once), thus allowing parallel edges between pairs of vertices. For 
simplicity we will imprecisely call G a graph and £ the edge set, even if G contains self-loops 
or multiple edges. The DaC model on a general (finite or infinite) graph G with vertex set V 
and edge set £, with parameters p, r G [0, 1], is defined by the following two-step procedure: 

• First step: Bernoulli bond percolation. We independently declare each edge in £ to be 
open with probability p, and closed with probability 1—p. We can identify a bond 
percolation configuration with an element rj £ {0, 1} £ : for all e G £, we define 77(e) = 1 if 
e is open, and 77(e) = if e is closed. 

• Second step: Bernoulli site percolation on the resulting cluster set. Given 77 £ {0,1} £ , we 
call p-clusters or bond clusters the connected components in the graph with vertex set V 
and edge set {e G £ : 77(e) = 1}. Note that the set of p-clusters of 77 gives a partition of V. 
For each p-cluster C, we assign the same color to all the vertices in C. The chosen color is 
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black with probability r and white with probability 1 — r, and this choice is independent 
for different p-clusters. 



These two steps yield a site percolation configuration £ G {0, 1} V by defining, for each v G V, 
= 1 if v is black, and £(v) = if v is white. The connected components (via the edge set £) 
in £ of the same color are called (black or white) r-clusters. The resulting measure on {0, 1} V 
is denoted by pi9 r . 

Let E^ C {0, 1} V denote the event that there exists an infinite black r-cluster. By standard 
arguments (see Proposition 2.5 in [S]), for each p G [0, 1], there exists a critical coloring value 
r£(p) G [0, 1] such that 



G r-pb 



6x1 = if r < rf(p), 
'• N 'l « if r - r< : {,>). 



This value should not be confused with the critical edge parameter pj? G [0,1], which is 
defined as follows: the probability that there exists an infinite bond cluster is for all p < p^, 
and positive for all p > p^ . It immediately follows from ergodicity that the latter probability 
is, in fact, 1 for all p > pjf, whence rf(p) = for all such p. Note, however, that although 
/i^^-E*-,) G {0, 1} for all p < p^ , there exist graphs (such as the square lattice, defined below) 
with n~ (E^q) G (0, 1) for some r > rf(p). 

Our main goal in this paper is to understand how the critical coloring parameter r^ depends 
on the edge parameter p. Since the addition or removal of self-loops obviously does not affect 
the value of rf(p), we will assume that all the graphs G that we consider in this paper are 
without self-loops. On the other hand, G is allowed to contain multiple edges. 

Our first result gives bounds on rf(p) in terms of (0), which is simply the critical value 
for Bernoulli site percolation on G. By the degree of a vertex v, we mean the number of edges 
incident on v (counted with multiplicity). 

Proposition 1. For any graph G with maximal degree A, it holds for all p G [0, 1) that 

l-rf(0) C/ . jf(0) 



(l_p)A - csr,- (i_ p) A- 

Proposition [1] will be a simple corollary of the much more general Theorem II T\ to which we 
will return in Section [5J 

An important special case arises when we take G to be the square lattice, whose vertex set 
is Z 2 , and edge set £ 2 is given by the edges between elements of Z 2 at Euclidean distance 1 
from each other. With an abuse of notation, we denote this graph by Z 2 . It has been shown 
in [H [8] that for all p > 1/2, we have rf (p) = 0, that rJ 2 (l/2) = 1, and that for p < 1/2, 
r ? (p) ^ [l/2j !)• In the interval p G [0, 1/2), the exact value of rf (p) is unknown. 

Recall that (0) is the critical value for Bernoulli site percolation on the square lattice. It is 
generally believed that rf (0) 0.593, but the best rigorous lower bound known to date, given 
in [5], is 0.556, and the best rigorous upper bound is 0.679, see [TTj. The previous proposition, 
combined with these bounds on r^ 2 (0), gives rigorous bounds on r^ 2 (p). 



Corollary 2. For all p G [0, 1) 



0.444 -2. . 0.679 
1 < r v) < 

(i-pY ~ rc [p) - (i-pY 
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Note that these bounds are new only for small values of p. In particular, the upper bound 
is nontrivial only for p < 1 — (0.679) 1 / 4 ~ 0.092. Since (jp) > 1/2 is known, the lower bound 
is interesting only where it implies rj? (p) > 1/2. This is the case for p < 1 — (0.888) 1 / 4 ~ 0.029. 

Our other main results are concerned with the monotonicity and continuity properties of 
the critical value function. Proposition 2.7 in [8j (or the more general Remark [16] in this 
paper) implies that when G is a tree, r^(p) is non- increasing on [0, 1], and continuous on [0, 1). 
Continuity in 1 fails, for example, for G = Z {i.e., the infinite regular tree of degree 2), for which 
r c(P) = 1 for all p < 1, but r^(l) = 0. A less trivial example for discontinuous rf is G = Z 2 , 
where the above mentioned results in [U [8] imply that the function 2 (p) has a discontinuity at 
p = pf =1/2, and also that monotonicity fails since rf (1/2) > rjf (p) for all p ^ 1/2. Taking 
p = 1/2 is "cheating" in that 1/2 is precisely the critical value for Bernoulli bond percolation on 
Z 2 , and the structure of bond clusters at 1/2 is fundamentally different from those for p < 1/2. 
The relevant questions, therefore, are whether monotonicity and continuity hold for p < p^f. 

The problem of continuity is solved at p = when the graph G has bounded degree: 

Proposition 3. For any graph G with bounded degree, r~(p) is continuous in p at 0. 

The above result is a simple consequence of Proposition [H which in turn relies on a stochastic 
domination argument. Since bond clusters can be arbitrarily large for p > 0, we cannot extend 
this argument to the interval [0,p^). In fact, we will show that Proposition [3] is sharp both in 
that there exists a graph of unbounded degree whose critical value function is discontinuous at 
< Pc , and also that continuity may fail for bounded degree graphs at any other place than 
(but, with the exception of 1, still in the subcritical regime). 

Proposition 4. There exists a graph G with pf > such that r^ is discontinuous at 0. 

Theorem 5. For all po G (0,1), there exists a graph G of bounded degree such that rf(p) is 
discontinuous at po < p~? . 

The proofs of these and all further results in this section are given in Section [H 
The construction that we will use to prove Theorem O is admittedly somewhat artificial, and 
the question remains whether there exist more regular, such as transitive or quasi-transitive, 
graphs with a discontinuity of r^ below p^. For the definition of (quasi-)transitivity, see, for 
example, Definition 1.2 in [9]. 

Open question 6. Is r^(p) continuous in p on [0,p^) for every quasi-transitive graph G? 

Our next result shows that not only near-trivial reasons (such as the essentially different 
structure of bond clusters in different phases) can cause non-monotonicity of the critical value 
function, moreover, that such a phenomenon may occur on a "nice" graph as well. 

Proposition 7. There exists a quasi-transitive graph G such that r^ is not monotone on the 
interval [0,p^). 

The graph that we define in the proof of Proposition [7] is quasi-transitive, but not transitive. 
In fact, all transitive graphs whose critical value functions were studied so far have been proved 
or conjectured to possess the monotonicity property. 

Open question 8. Is r^(p) a monotone function of p for p € [0,p^) for all transitive graphs 
G? 

Note that strict monotonicity cannot be expected for all transitive graphs as it fails for the 
triangular lattice T, where rj(p) = 1/2 for all p < pj (see [1]). Our feeling is, however, that 
Open question [8] can be answered in the affirmative. 
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Now we turn our attention from general graphs to the square lattice, whose critical value 
function, as mentioned above, is not monotone on [0,1] and has a discontinuity at 1/2. Our 
next result, Theorem below, implies that this is the only discontinuity point of rf 2 (p). The 
proof of this theorem will be given in Section 14.31 

Theorem 9. The critical coloring value rf 2 (p) is a continuous function of p on the interval 

pe [0,1/2). 

Unfortunately, we have not yet been able to prove (or disprove) monotonicity of rj 2 (p) for 
subcritical p. However, our numerical experiments suggest (see Conjecture I2"2"j) that (p) is 
strictly decreasing on the interval [0, 1/2). 

A brief outline of the paper is as follows. We set the notation and collect a few results 
from the literature in Section [TJ In Section [2J we stochastically compare r with Bernoulli site 
percolation in Theorem [TTT and show how this result implies Proposition [TJ We determine the 
critical value function for a class of tree-like graphs in Section and use these results for proving 
several of our main results in Sections 14.11 and 14.21 We then prove Theorem [9] in Section 14.31 
by a finite-size argument. Finally, Section O deals with the asymptotic behavior of r^(p) as p 
tends to pf. 

1 Definitions and notation 

We consider a graph G = (V, £) where the set of vertices V = {vo,vi,V2, ■ ■ ■} is countable. We 
define a total order "<" on V by saying that Vi < Vj if and only if % < j. In this way, for any 
subset V C V, we can uniquely define min(U) G V as the minimal vertex in V with respect 
to the relation "<". For a set S, we denote {0, 1} S by fig. We call the elements of Og bond 
configurations, and the elements of fly site configurations. As defined in the Introduction, in 
a bond configuration 77, an edge e G £ is called open if 77(e) = 1, and closed otherwise; in a 
site configuration £, a vertex v G V is called black if £(e) = 1, and white otherwise. Finally, 
for r] G fig and v G V, we define the bond cluster Q v (r/) of v as the maximal connected induced 
subgraph containing v of the graph with vertex set V and edge set {e G £ : 77(e) = 1}, and 
denote the vertex set of 6^(77) by C v (rj). 

For a G [0, 1] and a set 5, we define as the probability measure on fig that assigns to 
each s G S value 1 with probability a and with probability 1 — a, independently for different 
elements of S. We define a function 

<3? : fig x fl v — > J7 £ x fiy, 
(77, «) ^ (?7,£), 

where £(u) = K(min(C v (77))). For p, r G [0, 1], we define P^ r to be the image measure of <g) i/; 
by the function $, and denote by /Li^ r the marginal of P^ r on fly. Note that this definition of 
fip r is consistent with the one in the Introduction. 

Finally, we give a few definitions and results that are necessary for the analysis of the DaC 
model on the graph Z 2 . The matching graph Z 2 of the square lattice is the graph with vertex set 
Z 2 and edge set £ 2 = {(v, w) : v = (vi,V2),w = (wi, u>2) G Z 2 , max(|?;i — w\\, 1^2 — ^21) = !}• In 
the same manner as in the Introduction, we define, for a color configuration £ G {0, 1} Z , (black 
or white) *-clusters as connected components (via the edge set £ 2 ) in £ of the same color. We 
denote by &*(p, r) the P? r -probability that the origin is contained in an infinite black *-cluster, 
and define 

r*(p) = sup{r : @*{p,r) = 0} 



4 



for all p G [0, 1]. (Note that this value may differ from r c *{p)\) The main result in [1] is that 
for all p G [0, 1/2), the critical values rj 2 (p) and r*{p) satisfy the duality relation 

rf (p) + r*(p) = l. (1) 

We will also use the famous exponential decay result for subcritical Bernoulli bond per- 
colation on 1? . Let denote the origin in Z 2 , and for each n G N = {1,2,...}, let us 
define S n = {v G Z 2 : dist(v,0) = n} (where disi denotes graph distance), and the event 
M n = {rj G r? £ 2 : there is a path of open edges in rj from to S n }. Then we have the following 
result: 

Theorem 10 (|13[[T]). For p < 1/2, there exists ip(p) > such that for all n G N ; we have that 

vf{M n )<e~ n ^\ 

2 Stochastic domination 

In this section, we will prove Proposition[T]via a stochastic comparison between the DaC measure 
and Bernoulli site percolation. Before stating the corresponding result, however, let us recall 
the concept of stochastic domination. 

We define a natural partial order on by saying that £' > £ for £, £' G f2y if, for all 
» £ V, ^'(tO > £( v )- A random variable / : ily — >• M is called increasing if £' > £ implies that 
/(£') > /(£)> an d an event E C fiy is increasing if its indicator random variable is increasing. For 
probability measures fx, ll' on fly, we say that fx' is stochastically larger than fx (or, equivalently, 
that [x is stochastically smaller than //, denoted by ll < s t //) if, for all bounded increasing 
random variables / : fly —> R, we have that 

/ /(e) dfi'io > [ f(o dLi(o- 

By Strassen's theorem [16], this is equivalent to the existence of an appropriate coupling of 
the measures ll' and fx; that is, the existence of a probability measure Q on fly x fly such 
that the marginals of Q on the first and second coordinates are ll' and ll respectively, and 

Q({(e / ,e)G^vxf) v :e'>e}) = i. 

Theorem 11. For any graph G = (V, £) whose maximal degree is A, at arbitrary values of the 
parameters p,r G [0, 1], 

V G V 

U r(l-p) A - st ^P,r -st ^l-(l-r)(l-p) A - 

Before turning to the proof, we show how Theorem [11] implies Proposition [1] 

Proof of Proposition [1] It follows from Theorem [11] and the definition of stochastic dom- 
ination that for the increasing event E^ (which was defined in the Introduction), we have 
Hp ^E 1 ^) > whenever r(l — p) A > r~ (0), which implies that r~(p) < r^(0)/(l — p) A . The 
derivation of the lower bound for r^f(p) is analogous. □ 

Now we give the proof of Theorem [11] which bears some resemblance with the proof of 
Theorem 2.3 in [3]. 
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Proof of Theorem II 11 Fix G = (V, £) with maximal degree A, and parameter values 
p,r£ [0, 1]. We will use the relation "<" and the minimum of a vertex set with respect to this 
relation as defined in Section [TJ In what follows, we will define several random variables; we 
will denote the joint distribution of all these variables by P. 

First, we define a collection {jf X) y ■ x, y £ V, e = (x, y) £ £) of i.i.d. Bernoulli(p) random 
variables (i.e., they take value 1 with probability p, and otherwise); one may imagine having 
each edge e G £ replaced by two directed edges, and the random variables represent which of 
these edges are open. We define also a set (k x : x £ V) of Bernoulli(r) random variables. Given 
a realization of ijf x ^ '■ x, y £ V, e = (x, y) G £) and (k x : x £ V), we will define an x fig-valued 
random configuration with distribution P^ r , by the following algorithm. 

1. Let v = min{x £ V : no £-value has been assigned yet to x by this algorithm}. (Note 
that v and V,Vi,Hi (i £ N), defined below, are running variables, i.e., their values will be 
redefined in the course of the algorithm.) 

2. We explore the "directed open cluster" V of v iteratively, as follows. Define vq = v. 
Given vo, v\, . . . , V{ for some integer i > 0, set 77(e) = r\% i w for every edge e = (vi,w) G £ 
incident to Vi such that no 77-value has been assigned yet to e by the algorithm, and write 
Hi + i = {w £ V \ {vq, vi, . . . , V{} : w can be reached from any of Vq, vi, ■ ■ ■ , Vi by using 
only those edges e £ £ such that 77(e) = 1 has been assigned to e by this algorithm}. If 
Hi + \ ^ 0, then we define Vi + \ = min(i7j + i), and continue exploring the directed open 
cluster of v; otherwise, we define V = {vq, v\, . . . , V{}, and move to step 3. 

3. Define £(w) = k v for all w £ V, and return to step 1. 

It is immediately clear that the above algorithm eventually assigns a £-value to each vertex. 
Note also that a vertex v can receive a £-value only after all edges incident to v have already 
been assigned an 77-value, which shows that the algorithm eventually determines the full edge 
configuration as well. It is easy to convince oneself that (77, £) obtained this way indeed has the 
desired distribution. 

Now, for each v £ V, we define Z(v) = 1 if k v = 1 and 77^ v = for all edges e = (v, w) £ £ 
incident on v (i.e., all directed edges towards v are closed), and Z(v) = otherwise. Note 
that every vertex with Z(v) = 1 has £(v) = 1 as well, whence the distribution of £ (i.e., fJ,9 r ) 
stochastically dominates the distribution of Z (as witnessed by the coupling P). 

Since the distribution of Z is a product measure on f^v with parameter r(l — p) d ^ at v, 
where d{v) < A is the degree of v, we conclude that /i^ r stochastically dominates the product 
measure on fly with parameter r(l — p) A , which gives the desired stochastic lower bound. 
The upper bound can be proved analogously; alternatively, it follows from the lower bound by 
exchanging the roles of black and white. □ 

Remark 12. The same stochastic bounds hold for any DaC(q) measure [jS^^ (which is obtained 
by replacing Bernoulli bond percolation in the definition of the DaC model with a so-called 
random- cluster measure with parameters p and q) with q > 1, while one needs to replace 
p in both bounds by p/(p+ (1 — p)q) in case of q < 1 (assuming that the DaC(q) model on G 
exists in the latter case). Of course, these stochastic bounds also imply deterministic bounds for 
the critical coloring value in terms of (0). 

Sketch of proof. This is just a straightforward generalization of the above proof. We define 
an fiv x fig-valued random configuration (77, £) by an algorithm as above, with two differences. 
The first one is that here we define (77^^ : x, y £ V, e = (x, y) £ £) to be i.i.d. random variables 
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with uniform distribution on the interval [0,1]. Secondly, in step 2, we fix rj for those edges e 
incident to V{ that are without an //-value one at a time, taking 77(e) = for e = (vi, w) G £ 
if and only if 77^. w < (j), where cp is the conditional probability of e being closed under the 
conditional distribution ^ given the state of those edges whose 77-value has been determined 
before; otherwise, we take 77(e) = 1. This indeed gives the desired distribution: since the 
algorithm assigns ^-values only after the whole bond cluster of the vertices in question has been 
determined, the ^-values do not affect the further edge distribution. 

Since the ^^-probability of an edge e being closed given the state of every other edge is 
either 1 — p or 1 — p/(p + (1 — p)q) (see, e.g., [2]), it follows that 4> > b with 6=1— p i£ q > 1 
and 6 = 1— p/(p + (1 — p)q) if q < 1. Now define, for each vertex v, Z(v) = 1 if k v = 1 and 
77^ v < b for all edges e = (v, w) G £ incident on v, and Z(v) = otherwise. Upon noting that 
the distribution of Z is a product measure on with parameter rb d ^ at v and that 77^, v < b 
implies that the corresponding directed edge from w to v is closed, the rest of the proof goes as 
above. □ 

3 The critical value functions of tree-like graphs 

In this section, we will study the critical value functions of graphs that are constructed by 
replacing edges of an infinite tree by a sequence of finite graphs. We will then use several such 
constructions in the proofs of our main results in Section 01 

Let us fix an arbitrary sequence D n = (V n ,£ n ) of finite connected graphs and, for every 
n G N, two distinct vertices a n ,b n £ V n . Let I^ = (V^-E^) denote the (infinite) regular tree 
of degree 3, and fix an arbitrary vertex p G V3. Then, for each edge e G E%, we denote the 
end-vertex of e which is closer to p by /(e), and the other end-vertex by s(e). Let Yd = (V,E) 
be the graph obtained by replacing every edge e of I^ between levels n — 1 and n (i.e., such 
that dist(s(e), p) = n) by a copy D e of D n , with a n and b ri replacing respectively /(e) and s(e). 
Each vertex v G V3 is replaced by a new vertex in V, which we denote by v. It is well known 
that p^ 3 = r^ 3 (0) = 1/2. Using this fact and the tree-like structure of Fd, we will be able to 
determine bounds for p\ D and r\ D (j>). 

First, we define h D "(p) = v^ n {a n and b n are in the same bond cluster), and prove the fol- 
lowing, intuitively clear, lemma. 

Lemma 13. For any p G [0, 1], the following implications hold: 

a) if lim sup^^ h Dn (p) < 1/2, thenp<p T c D ; 

b) iflimmf n ^ OD h Dn (p) > 1/2, then p > p T c D . 

Proof. We couple Bernoulli bond percolation with parameter p on with inhomogeneous 
Bernoulli bond percolation with parameters h Dn (p) on T3, as follows. Let 77 be a random variable 
with law u^, and define, for each edge e G £3, W(e) = 1 if /(e) and s(e) are connected by a 
path consisting of edges that are open in 77, and W(e) = otherwise. The tree-like structure 
of Yd implies that W(e) depends only on the state of the edges in D e , and it is clear that if 
dist(s(e), p) = n, then W(e) = 1 with probability h Dn (p). 

It is easy to verify that there exists an infinite open self-avoiding path on Yd from p in the 
configuration 77 if and only if there exists an infinite open self-avoiding path on T3 from p in the 
configuration W. Now, if we assume lim sup n _> 00 h Dn (p) < 1/2, then there exists t < 1/2 and 
N G N such that for all n> N, h Dn (p) < t. Therefore, the distribution of the restriction of W 
on L = {e G E3 : dist(s(e), p) > N} is stochastically dominated by the projection of vf A on L. 
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This implies that, a.s., there exists no infinite self-avoiding path in W, whence p < p\ D by the 
observation at the beginning of this paragraph. The proof of b) is analogous. □ 

We now turn to the DaC model on To. Recall that for a vertex v, C v denotes the vertex 
set of the bond cluster of v. Let E an ^ n C Og n x J7v n denote the event that a n and b n are 
in the same bond cluster, or a n and b n lie in two different bond clusters, but there exists a 
vertex v at distance 1 from C an which is connected to b n by a black path (which also includes 
that £(v) = £,{b n ) = 1). We note that this is the same as saying that C an is pivotal for the 
event that there is a black path between a n and b n , i.e., that such a path exists if and only 
if C an is black. It is important to note that E an ^ n is independent of the color of a n . Define 
f Dn (p,r) = F®™(E antbn ), and note also that, for r > 0, f Dn {p,r) = P^, 1 (there is a black path 
from a n to b n | £(a„) = 1). 

Lemma 14. For any p,r £ [0, 1], we have the following: 

a) if limswp n _^. O0 f Dn (p, r) < 1/2, then r < r^ D (p); 

b) if lim inf n ^oo f Dn (p, r) > 1/2, then r > r^ D (p). 



Proof. We couple here the DaC model on T/j with inhomogeneous Bernoulli site percolation 
on IV For each v € V3 \ {p}, there is a unique edge e £ E3 such that v = s(e). Here we denote 
D e (i.e., the subgraph of replacing the edge e) by Dy, and the analogous event of E ari ^ n for 
the graph D v by E%. Let (77, £) with values in fig x Qy be a random variable with law 
We define a random variable X with values in Oy 3 , as follows: 



X(v) 



£(p) ifv = p, 

1 if the event E% is realized by the restriction of (77, £) to Di 
otherwise. 



As noted after the proof of Lemma [TBI if u = f((u, v)), the event E^ is independent of the 
color of u, whence (Ey) vG y 3 \^ p j are independent. Therefore, as X(p) = 1 with probability r, 
and X(v) = 1 is realized with probability f Dn (p,r) for v € V3 with dist(v,p) = n for some 
n £ N, X is inhomogeneous Bernoulli site percolation on 

Our reason for defining X is the following property: it holds for all v E V3 \ {p} that 

p -R- v if and only if p «->• v, (2) 

z 

where x y denotes that x and y are in the same black cluster in the configuration Z. Indeed, 

assuming p •<-)■ v, there exists a path p = xq, x\, ■ ■ ■ , x\. = v in such that, for all < i < k, 

Xi A x£±.i holds. This implies that £(/3) = 1 and that all the events (i?x- i )o<j<fe occur, whence 

X(xi) = 1 for i = 0, . . . , k, so p •<-)■ v is realized. The proof of the other implication is similar. 
It follows in particular from @ that p lies in an infinite black cluster in the configuration £ if 
and only if p lies in an infinite black cluster in the configuration X. 

Lemma [T^l presents two scenarios when it is easy to determine (via a stochastic comparison) 
whether the latter event has positive probability. Assume, for example, that liminf n _ !>00 f Dn (p, r) > 
1/2, whence there exists t > 1/2 and N £ N such that for all n > N, f Dn (p, r) > t. In this case, 
the distribution of the restriction of X on K = {v £ V3 : dist(v, p) > A^} is stochastically larger 
than the projection of vf 3 on K. Let us further assume that r > 0. In that case, X(p) = 1 
with positive probability, and f Dn (p, r) > for every n £ N. Therefore, under the assumptions 
liminfn^oo f Dn (p,r) > 1/2 and r > 0, p is in an infinite black cluster in X (and, hence, p is in 
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an infinite black cluster in £) with positive probability, which can only happen if r > r^ D (p). On 
the other hand, if liminf n _ ) . 00 f Dn (p, 0) > 1/2, then it is clear that lim uhv^oo f Dn (p, r) > 1/2 
(whence r > r^ D (p)) for all r > 0, which implies that r^ D (p) = 0. The proof of part a) is 
similar. □ 



When D n = D for some finite connected graph D = (V,E), the previous lemmas determine 
the critical p- and r-values, as follows. Since D is finite and connected, h D (p) is a polynomial 
in p with h D (0) = and h D (l) = 1, and it is strictly increasing on [0,1]. We can therefore 
define po £ (0,1) as the unique value such that h D (po) = 1/2. Also, for fixed p G [0,1], 
f D (p,0) = h D (p) and f D (p,l) = 1, and it is also clear that, for any fixed p < 1, f D (p,.) is 
strictly increasing and continuous on [0, 1]. Hence, for p < po, we can define rn{p) G (0, 1) to 
be the unique value such that f D (p,rD(p)) = 1/2. For p > pjj, we define rn(p) = 0, and note 
that f D {p,r) > 1/2 for all r > 0. 

Lemma 15. If D n = D, then p T c D = po, and for all p G [0, 1], r^ D (p) = rr>(p)- 

Proof. These statements immediately follow from Lemmas ll3l and ll4l upon noting that lim ra _ il00 
h D (p), and similarly with f D ". □ 

Remark 16. A trivial extension of the above proofs gives that ifT^ is replaced in the definition 
ofTjj by an arbitrary tree V, then analogous statements to those in Lemmas \13H15\ hold with 
1/2 replaced everywhere by l/br(T), where br(T) is the branching number ofT (see ]12^ for the 
definition). Furthermore, ifTz> is constructed by taking an arbitrary enumeration {e±,C2, ■ ■ ■} 
of the edge set of T and replacing e« by D{ for every ieN, the same implications still hold. 

4 Proofs of the main results 

We will first apply the machinery developed in Section [3] for proving some of our main results 
in Sections 14. II and 14.21 In these sections, we will freely use the notations introduced in Section 
[3J . We then turn to the proof of Theorem [9] in Section 14.31 

4.1 Non-monotonicity 

The results in Section [3] enable us to prove that (a small modification of) the construction 
considered by Haggstrom in the proof of Theorem 2.9 in [8] is a graph whose critical coloring 
value is non-monotone in the subcritical phase. 

Proof of Proposition 0. We will apply Lemma [T5l with D n = D = D k , which is defined, 
for k £ N, to be the complete bipartite graph with the vertex set partitioned into {^1,^2} 
and {a,b, v\,V2, ■ ■ ■ ,Vk} (see Figure [TJ. We call ex,e[ and e2,e' 2 the edges incident to a and b 
respectively, and for i = 1, . . . , k, fi, f[ the edges incident to Vi. For each k £ N, we denote the 
graph constructed with D k by T D k. The resulting graphs, as is always the case when D n is the 
same for all n, are clearly quasi-transitive. 

We will show below that it holds for all k G N that 



p D k > 1/3, 

7-^(0) < 2/3, and 

r D „(l/3) < 2/3. 



(3) 
(4) 
(5) 
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Figure 1: The graph D k . 



Furthermore, there exists k £ N and po £ (0, 1/3) such that 

r Dk ( Po )> 2/3. (6) 

Proving ([3])-([6]) will finish the proof of Proposition [7] since these inequalities imply by LemmafTBI 
that, with k £ N as in (0), the quasi-transitive G = T D k is a graph with a non- monotone critical 
value function in the subcritical regime. 

Throughout this proof, we will omit superscripts in the notation when no confusion is 
possible. For the proof of ([3]), recall that hP is strictly increasing in p, and h (PD k ) = V^- 
Since 1 — h (p) is the z/ p -probability of a and b being in two different bond clusters, we have 
that 

1 - hP (1/3) > 1/1/3 ({ei and e'i are closed} U {e 2 and e' 2 are closed}). 

From this, we get that h Dk (1/3) < 25/81, which proves ([3]). 

To get we need to remember that for fixed p < Pi)k, f D (p,r) is strictly increasing in 
r, and f Dk (p,r D k(p)) = 1/2. One then easily computes that /(0, 2/3) = 16/27, whence (jl]) 
follows. 

Now, define A to be the event that at least one edge out of e\, e' 1; e 2 and e' 2 is open. Then 

/ Dfc (l/3,2/3) > Pl/ 3 ,2/3(^a,6 | A)F 1/3) 2/ 3 (A) 

> Vi/3,2/a(C b black | A) ■ 65/81, 

which gives that / Dfc (l/3,2/3) > 130/243, and implies ©. 

To prove ©, we consider to be the event that ei, e' l5 e 2 and e' 2 are all closed and that 
there exists i such that fi and f- are both open. One can easily compute that 

P p , r ( J B fc ) = (l-p) 4 (l-(l-p 2 ) fc ), 

which implies that we can choose po £ (0, 1/3) (small) and k £ N (large) such that ¥ po ^(Bi c ) > 
17/18. Then, 

f D \p , 2/3) = P po , r (£ a , fe | S fc )P P0 , r (B fc ) + F P0 A E a,b I ^)(l-F POir (B fc )) 
< (2/3) 2 • 1 + 1- 1/18, 

whence inequality ([6]) follows with these choices, completing the proof. □ 
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4.2 Graphs with discontinuous critical value functions 

Proof of Proposition For n G N, let -D n be the graph depicted in Figure [H and let G be 

constructed with this sequence of graphs as described at the beginning of Section [3l 




Figure 2: The graph D 



It is elementary that lim n _ il00 hP n (p) = p, whence pj? = 1/2 follows from Lemma flUl thus 
p = is subcritical. Since lim n _ i , 00 f Dn (0,r) = r 2 , Lemma fbH gives that (0) = 1/V2. On the 
other hand, lim n _). 00 f Dn (p, r) = p + (1 — p)r for all p > 0, which implies by Lemma [bH that for 
P < 1/2, 

r , » 1/2 -p 

as p — > 0, so is indeed discontinuous at < □ 

In the rest of this section, for vertices v and w, we will write v w to denote that there 
exists a path of open edges between v and w. Our proof of Theorem [5] will be based on the 
following modification of Lemma 2.1 in |14j : 

Lemma 17. For all po,p\ G (0,1), f/iere exists a sequence G n = (V n ,E n ) of graphs and 
XniVn G V n of vertices (n G N) such that 

1. Vp ( "(x n O y n ) > pi /or a// n; 

2. lim^oo Vp n (x n H!/„) = for all p < p , and 

3. there exists A = A(po>Pi) < o° such that, for all n, G n has maximal degree A. 



Proof. There are several ways of constructing such a sequence. Since the concrete choice of 
G n will be irrelevant for the rest of this paper, we will just show how to modify the proof of 
Lemma 2.1 in [T3], using the terminology of that proof. 

Choose £2 > and c > such that (1 — £2) 2 (1 — e _c / 3 + 1 ) > Pl) fix e± = £2/3, and let m be 
so large that Vp (0 is in an infinite bond cluster) > a/1 — e\. Fixing A n = {v{ : 1 < i < c/pq}, 
it is easy to see that Vn 3H n : Vt> G A n ^ n (0 f> v) > 1 — E\. Defining G n with these H n as in 
[LI] and taking z n and y n as in if n and the copy of H n respectively, the computations in the 
proof of [2] give the statement of Lemma [TTJ D 

Lemma [T7] provides a sequence of bounded degree graphs that exhibit sharp threshold- type 
behavior at pq. We will use such a sequence as a building block to obtain discontinuity at po i n 
the critical value function in the DaC model. The basic idea behind the following construction 
was suggested to us by Jeff Steif. 



Proof of Theorem O Fix Po G (0, 1) and the smallest k G N such that 

(i±^)'<v, 
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As ((1 + 2p )/3) fc_1 > 1/2, we may fix 




and also G n = (V n ,E n ),x n ,y n (n G N) for po and p\ as in Lemma [T71 Let P be a graph with 
vertex set {vq, v±, . . . , Vk+i} and edge set {(vi,Vj) : \i — j\ = 1} (i.e., P is a path with k edges). 
We define, for each n G N, D n by replacing the edge between Vk and t>fc + i by G n , with x n 
replacing vj. and y ra replacing Ufc+i. Let a n and 6 n be the two vertices of D n that correspond 
to the former vo and Vk+i, respectively, and let G be the graph Yd with the sequence D n as 
defined in Section [3j 

Since h Dn (p) = p k v^ n (x n «-»■ y n ) < p fc for all n, we have limsup /i Dn ((l + 2p )/3) < 1/2, 
whence Lemma [13] implies that po < (1 + 2po)/3 < p^ . Define 

Jo if k = 1; 

C= \(^^H iffc>l, 

I l-Po ' 

and note that c < 1/2 by the choice of p\. We will show below that r^(po) < c while r^(p) > 1/2 
for all p < Po, which implies discontinuity of r£* at po < p^ , finishing the proof. 
Fix an arbitrary r G (c, 1/2), and note that for all p G [0, 1], 

f D -(p,r) = (p+(l-p)r) k - l (vf(x n o y n ) + (1 - vf(x n y n ))a), 

with a being the conditional -probability of x n and y n being connected by a black path, 
given that they are not in the same bond cluster and x n is black. It is immediately clear that 
a < r, whence we have for any fixed p < po that lim sup n _ >00 f Dn (p, r) < lirrin^oo (x n o 
Vn) + r < 1/2 by r < 1/2, which implies by Lemma [bH that rf(p) < r. On the other hand, 
liminf^oo f Dn (p ,r) > (p + (1 - Po)^' 1 liminf n _ ) . 00 v^ n [x n O y n ) > 1/2 by r > c, whence 
Lemma [TH gives that r > r^(po). These computations go through for any r G (c, 1/2), thus we 
indeed have rf{p) > 1/2 for all p < po and rf(po) < c. □ 

4.3 Continuity of rjf 2 (p) on the interval [0,1/2) 

In this section, we will prove Theorem [9j Our first task is to prove a technical result valid on 
more general graphs stating that the probability of any event A whose occurrence depends on a 
finite set of £- variables is a continuous function of p for p < p^ . The proof relies on the fact that 
although the color of a vertex v may be influenced by edges arbitrarily far away, if p < p^?, the 
corresponding probability decreases to in the limit as we move away from v. Therefore, the 
occurrence of the event A depends essentially on a finite number of rj- and ft- variables, whence 
its probability can be approximated up to an arbitrarily small error by a polynomial in p and 
r. 

Once we have proved Proposition [18] below, which is valid on general graphs, we will apply 
it on Z 2 to certain "box-crossing events," and appeal to results in [3j to deduce the continuity 
of (p). 



Proposition 18. For every site percolation event A C {0, 1} V depending on the color of finitely 
many vertices, Hp r (A) is a continuous function of (p,r) on the set [0,pj?) x [0, 1]. 
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Proof. In this proof, when \i is a measure on a set S, X is a random variable with law [i and 
F : S — > R is a bounded measurable function, we write abusively fi[F(X)] for the expectation 
of F(X). We show a slightly more general result: for any k > 1, x = (x%, . . . , x^) G V fc 
and / : {0, l} fc — > R bounded and measurable, [i^ T [/(£(xi), . . . ,£(xk))] is continuous in (p,r) 
on the product [0,p^) x [0,1]. Proposition [18] will follow by choosing an appropriate family 
{x\, . . . such that the states of the X{ suffices to determine whether A occurs, and take / 
to be the indicator function of A. 

To show the previous affirmation, we condition on the vector m x {rj) = (min C Xl (rf), ... , min C Xk 
which takes values in the finite set 

V = . . . , v k ) G V k : Vivi < max{xi, . . . , 
and we use the definition of P^ r as an image measure. By definition, 
Mg r [/(^l),---,fe))] 

= E ® ^ • • • ' K ( u fc))l{ m a3 = «}] [{ m ^ = «}] 



Note that z/ T v [/(k(ui), • • • , «(i>fc))] is a polynomial in r, so to conclude the proof we only need to 
prove that for any fixed x and v, Up ({ m(x) = v }) depends continuously on p on the interval 
[0,p c G ). 

For n > 1, write = {|C a;i | < n, . . . , |C Xft | < n}. It is easy to verify that the event 
{m x = v} n F n depends on the state of finitely many edges. Hence, Up [{m x = v} n F n ] is 
a polynomial function of p. 

Fix po < Pc ■ F° r an P — POi 

< rf [{m(x) = v}} - rf [{m x = v} n F n ] < rf[F^} 

< 4 \K\ 

where lim v% [FS] = 0, since po < p*f . So, u„ [mix) = v] is a uniform limit of polynomials on 

any interval [0,poL Po < Pc> which implies the desired continuity. □ 

Remark 19. In the proof we can see that, for fixed p < p^ , fj,p >r (A) is a polynomial in r. 

Remark 20. If G is a graph with uniqueness of the infinite bond cluster in the supercritical 
regime, then it is possible to verify that v£ [{m(x) = v}] is continuous inp on the whole interval 
[0, 1]. In this case, the continuity given by the Proposition\T8l can be extended to the whole square 
[0,1] 2 . 

Proof of Theorem[9]. In order to simplify our notations, we write P Pi r, v p , respectively, r c (p), 
for fp 2 r , rf 2 and rf 2 (p). Fix po E (0, 1/2) and e > arbitrarily. We will show that there exists 
<5 = <5(pcb £ ) > such that for all p £ (po — d~,po + 5), 

r c {p) > r c (p ) - e, (7) 

and 

r c (p) < r c (p ) + e. (8) 
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Note that by equation ([T]), for all small enough choices of 5 > (such that < po db S < 1/2), 
([7]) is equivalent to 

r* c (p)<r*(p )+e. (9) 

Below we will show how to find 5\ > such that we have ((8j) for all p £ (po — Si , po + 6\) . One 
may then completely analogously find S2 > such that Q holds for all p 6 (po — ^2>Po + ^2)) 
and take <5 = min^i, £2). 

Fix r = r c (po) + £, and define the event V n = {(£,77) £ Slj? x Og 2 : there exists a vertical 
crossing of [0,n] x [0, 3n] that is black in £}. By "vertical crossing," we mean a self-avoiding 
path of vertices in [0, ra] x [0, 3n] with one endpoint in [0, n] x {0}, and one in [0, n] x {3n}. 
Recall also the definition of M n in Theorem 1101 By Lemma 2.10 in [3], there exists a constant 
7 > such that the following implication holds for any p, a £ [0, 1] and L £ N: if both 
(3L + 1)(L + l)z/p(M|£/ 3 j ) < 7 and P Pl a(V£,) > 1 — 7 are satisfied, then a > r c (p). (As usual, 
[x\ for x > denotes the largest integer m such that m < x.) Fix such a 7. 

By Theorem I10| there exists N £ N such that 

(3n+l)(n + l)^ (M Ln/3j ) < 7 

for all n> N. On the other hand, since r > r c (po)j it follows from Lemma 2.11 in [4j that there 
exists L > N such that 

Pp ,r(VL)>l- 7 . 

Note that both (3L + 1)(L + 1)i/ p (M\l/ s \ ) and P Pi r(VL) are continuous in p at po- Indeed, the 
former is simply a polynomial in p, while the continuity of the latter follows from Proposition [THJ 
Therefore, there exists 61 > such that for all p £ (po — <5i>R) + <$i)> 

(3L + l)(L + l)i/ p (M LV3 j) < 7, 

andPp,^^) > I-7. 

By the choice of 7, this implies that r > r c (p) for all such p, which is precisely what we wanted 
to prove. 

Finding 62 > such that © holds for all p £ (po — 82, Po + ^2) is analogous: one only needs 
to substitute r c (po) by r*(po) and "crossing" by "*-crossing," and the exact same argument as 
above works. It follows that 5 = min(<5i, 62) > is a constant such that both ([8]) and ([9]) hold 
for all p £ (po — 5, po + 5), completing the proof of continuity on (0, 1/2). Right-continuity at 
may be proved analogously; alternatively, it follows from Proposition [3J 

Remark 21. It follows from Theorem^ and equation (Tip that r*(p) is also continuous in p on 
[0,1/2). 

5 Asymptotic behavior as p tends to 

We performed an extensive simulation of the DaC model on the square lattice 1? and on the 
hexagonal lattice H, for various values of p; our results will be published separately in a separate, 
more numerically oriented paper [3]. In all this section, only these two lattices will be considered, 
and the letter G will be used to denote either of them. We obtained the confidence intervals 
represented in Figure [3j 

Having looked at Figure El we conjecture the following concerning the behavior of r^(p) as 
a function of p: 
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0.1 0.2 0.3 0.4 



0.1 0.2 0.3 0.4 0.5 0.6 0.7 



Figure 3: Simulation results for different values of p < p c (left: on the square lattice; right: on 
the hexagonal lattice). The dashed line was obtained via a non-rigorous correction method. 



Conjecture 22. For G denoting either Z 2 or H, in the interval p E [0,p c (G)), r^(p) is a 
strictly decreasing function ofp, and lim p ^. Pc (c)_ r c (p) = 1/2. 

Since it is rigorously known that r^(0) > 1/2 and r^{p) > 1/2 for all p G [0,p c (G)), 
Conjecture 1221 would imply that r^{p) > 1/2 for all p < p c (G). This suggests that the DaC(l) 
model on Z 2 or H is qualitatively different from the DaC(l) model on the triangular lattice, 
where the critical value of r is 1/2 for all subcritical p (see Theorem 1.6 in [3]). However, 
limp^x^- r c (p) = 1/2 would mean that the difference disappears as p tends to p c (G). 

The fact that the difference should disappear was conjectured by one of the authors (VB) 
and Federico Camia, based on the following heuristic reasoning. Near p = p c (G), the structure 
of the random graph determined by the bond configuration (whose vertices correspond to the 
bond clusters, and there is an edge between two vertices if the corresponding bond clusters 
are adjacent in Z 2 ) is given by the geometry of "near-critical percolation clusters," which is 
expected to be universal for 2-dimensional planar graphs. This suggests that the critical r for 
p close to its critical value should not depend much on the original underlying lattice, and we 
expect the convergence of r c (p) to 1/2 to be universal and hold in the case of any 2-dimensional 
lattice. 

There is an additional, strange feature appearing in the case of the square lattice: r c (p) 
seems to be close to being an affine function of p on the interval [0, 1/2). This is not at all the 
same on the hexagonal lattice, and we have not found any interpretation of this observation, or 
of the special role Z 2 seems to play here. 

Open question 23. Is r^" (p) an affine function ofp for p < 1/2? 

Acknowledgment. We thank Jeff Steif for suggesting (a variant of) the graph that appears 
in the proof of Theorem [5j 
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